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Abstract 

Short geodesies are important in the study of the geometry and the spec- 
tra of Riemann surfaces. Bers' theorem gives a global bound on the length of 
the first 3g— 3 geodesies. We use the construction of Brooks and Makover of 
random Riemann surfaces to investigate the distribution of short (< log(g)) 
geodesies on a random Riemann surfaces. We calculate the expected value 
of the shortest geodesic, and show that if one orders prime non-intersecting 
geodesies by length 71 < 72 < • • • < 7i, • • • , then for fixed k, if one allows 
the genus to go to infinity, the length of 7^ is independent of the genus. 

1 Introduction 

A standard tool in the study of compact Riemann surfaces is the decomposition 
into "pairs of pants" (Y pieces). Given a surface of genus g > 2, there are 3g — 
3 simple closed geodesies which partition the surface into g — 1 such pieces. 
Bounds on the lengths of the geodesies in such partitions are extremely desirable. 
If the geodesies are 71, . . . ,735-3, their lengths give half of the Fenchel- 
Nielsen parameters which parametrize the 6(7 — 6 dimensional Teichmiiller space 
of compact surfaces of genus g. 

Bers (0101) proved that for every compact Riemann surface of genus g > 2 
there is a partition with 
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where L g is a constant depending only on g. The best possible such constant 
is called Ber's constant. A constructive argument due to Abikoff (|1]) gives an 
explicit bound for L g ; unfortunately this bound grows faster than exponentially in 
g. The best result known is 

Theorem 1.1 (|10|). Every compact Riemann Surface of genus g < 2 has parti- 
tion 71, . . . , 73 9 -3 satisfying 

Klk) < 4k log 8?r( ^~ 1} k = 1, . . . , 3g - 3 

The longest geodesic is bounded by 

73,-3 < 26(0-1) (1) 

One might hope that in fact the bound in ([TJ might be improved to a logarith- 
mic one, but this is impossible. The "hairy torus" gives a lower bound of 

Proposition 1.1. L g > y/6g — 2 for all g > 2. 

The length of the shortest geodesic is also of particular interest. There 
are examples of classes of surfaces, such as the conformal compactification of the 
principal modular surfaces, where it is known that (J9j|7| ) 

Z( T i) = O(log0). 

Recently Katz Schaps and Vishne ifTol show for Hurwitz surfaces and some 
other principal congruence subgroups of arbitrary arithmetic surfaces a similar 
behavior. But all these examples are rare in the sense that they not occur in all 
genera and when they do occur they there are a small number of such surfaces. 
In this paper, we study the Belyi surfaces . Such surfaces are dense in the set of 
Riemann surfaces (|2|). In contrast to the previous examples , we prove 

Theorem 1.2. Let S be a Belyi surface of genus g > 2, as g — » oo the length of 
the shortest simple closed geodesic on S, denoted syst(S), is bounded by 

2.809 < E(syst(S)) < 3.085. (2) 

When E is the expected value. In particular, E(syst(S)) is independent of g for 
surfaces with large genus. 
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We also get information about the lengths of some longer geodesies. We con- 
sider a set of prime non intersecting geodesies arranged by length 71 < 72 < 
• ■ • < 7i, • • • • We will show that they are all simple and therefore, it is always 
possible to complete this set and get a pair of pants decomposition ifTOl . We show 
that 

Theorem 1.3. For any fixed i, E{l{^i)) is independent of the genus g as g — > 00. 

This estimate in fact holds for roughly short geodesies out of the total 
3(7 — 3 which requires to get a pants decomposition of the surface. We provide 
estimates for the lengths of these geodesies and show that they are all bounded 
by C log g. This result is in sharp contrast to Buser's result for general Riemann 
surfaces where the shortest geodesies are ~ log g, the same magnitude as the 
longest geodesies covered by our estimate. 

We analyze the length of the geodesies using the method developed by Brooks 
and Makover in JE1. They use cubic graph to generate the Belyi surfaces and en- 
dow the set surfaces with probability measure inherit from random cubic graphs. 

In section El we will describe the basic construction of compact surfaces from 
cubic graph. We review the connection between the metric structure of the surface 
and the combinatorics of the graphs. Next, in section |3] we connect cycles in the 
graph with geodesies on a surface, and show that our methods allow us to consider 
cycles in the graph Y of length ~ log g (seefTOI). We will do so by investigating 
an interesting explicit example of random matrix multiplication and analyzing the 
different moments of the distribution of the entries. As it turns out, the entries of 
the matrix product are a Stern sequence which is an integer sequence of indepen- 
dent interest. 

Next, in section HI (Theorem 14.21) we use the above results to show that the 
length of the geodesic associated to such a cycle is growing linearly with the 
length of the cycle. A detailed analysis of the growth rate of the Stern sequence 
shows that for a given cycle length on the graph the length of geodesies on the 
surface is concentrating around the mean length. We produce the Stern sequence 
using random products of 2 x 2 matrices. 

These random products are similar to processes described by Viswanath ( %UY ) 
for random Fibonacci sequence, Lima and Rahibe ifTTll in the physics of disordered 
systems, and others. In our example, due to the properties of the Stern sequence, 
we get more detailed results on the random product then we could find in the 
literature for similar processes. This example may be of an independent interest 
in the study of random products of matrices, hence we include a more thorough 
discussion of our investigations than is necessary to prove our main results. 
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Finally, in section |3] we calculate the numerical bound to length of the shortest 
geodesic. 

This gives a somewhat surprising picture of such random Riemann surfaces. In 
([8 1) Brooks and Makover showed that the Cheeger constant of random Riemann 
surfaces is bounded from below. Hence the geodesies that we find can not discon- 
nect large pieces of the surface. The picture becomes even more complicated by 
the result dEHHl) that shows that there is an embedded ball on the surfaces with 
area of ~ | of the area of the surface. A "generic" Riemann surface therefore 
looks something like figure [T] 




Figure 1: A "generic" Riemann surface. 



2 Basic Construction 

We wish to construct Riemann surfaces from cubic (3 regular) graphs[8| . We 
will need to add additional combinatorial structure an orientation. We define ori- 
entation O on the graph, as a cyclic permutation of the edges emanating from a 
each vertex v. Each vertex in the graph will be identified with an ideal hyperbolic 
triangle T with vertices at 0, 1, and oo (see Figure El. We may think of the points 
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Figure 2: An ideal triangle. 



in P = {i,i + 1, ^} as "midpoints" of the corresponding sides of the triangle 
T. The solid lines in Figure|2]are geodesies joining the points in P with the point 
1+ 2 ' while the dotted lines are horocycles joining pairs in S. Now we can glue 
neighboring triangles subject to two conditions the first is that the "midpoints" 
of sides are glued together, and the second is that the orientation O at the vertex 
agree with the orientation of T. 

Such a gluing is uniquely determined given (T, O). The resulting surface is 
an open complete Riemann surface S°(T, O), and we define the compact surface 
S C (T, O) as its conformal compactification. While for any given graph of size 
n there are only 2 n such surfaces, the following theorem will allow us to model 
"generic" surfaces using surfaces generated from graphs. 

Theorem 2.1 (|2|). The set of surfaces constructed from cubic graphs is dense in 
the set of Riemann surfaces. 

To work with the set of surfaces generated from random cubic graphs with n 
vertices, we will need a probability measure for the set of such graphs. While the 
fact that for each n we have a finite set of such graphs, it is hard to work with ll2"2l 
[6 1 therefore we will use Bollobas' configuration model. 
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Put 6n balls in a hat; label the balls using the numbers 1,2,..., (2n), with 
three copies of each number. Then pick at random pairs of balls. We can define a 
graph by taking a set of 1, 2, ... , (2n) vertices, and connecting vertices V{ to Vj by 
an edge if a pair of two balls marked with i and j have been picked together. We 
will endow the set of oriented cubic graphs with a probability measure by picking 
a graph using the configuration model and then flipping an unbiased coin at each 
vertex to pick an orientation. Note, we will allow loops and double edges since 
they not do not interfere with the construction of the surface S° (T,0). We use the 
notation T n for the set of cubic graphs on 2n vertices with the above probability 
measure, and J 7 * for the set of oriented cubic graphs with the same probability 
measure. 

For the unoriented graphs Bollobas proved (@ 

Theorem 2.2. Let X, denote the number of closed paths in V of length i. Then the 
random variables Xj on T n are asymptotically independent Poisson distributions 
with means 



In order to study the lengths of simple closed geodesies on the compact sur- 
faces S C (T, O), we must understand the relation of the metric structure of S°(T, O) 
and S°(T, O). The following two theorems show that as n — > oo, almost all sur- 
faces S C (T, O) have a global metric structure arbitrarily close to S°(T, O). 

Theorem 2.3 (| 7]). We say that S°(T,0) has cusps > L if there is a set of 
disjoint horocycles, one around each cusp and each has length > L. Then for 
every e, there exists numbers L, r, and y such that, if the cusps of S° have length 
> L, outside the union of cusp neighborhoods U = Uf =1 /~ 1 (C y ) C S° of the 
cusps C\ and V = U^ =1 B r (pi) C S c , the metrics ds^ and ds 2 Q satisfy 

-. rds% < ds% < (1 + e)ds%. 

(1 + e) u ~ ° - v ; u 

The large cusp condition is necessary and enables us to compare the global 
metric of open and compact surfaces. Let Q be a property of 3-regular graphs 
with orientation, and denote by Prob n [Q] the probability that a pair (T, O) picked 
from JF* has property Q. 

Theorem 2.4 (|8|). For every L > 0, as n — > oo, 

Prob n [S° (r, O) has cusps of length > L ] — > 1. 



6 



3 Geodesies on S C (T, O). 



The discussion in the previous section shows that we can use the graph T to get 
information about the surface S C (T, O). To do this, we must begin with a cycle 
in T and its' associated geodesic on S°(T, O). Next, we track the geodesies as 
the cusps of S°(T, O) are closed to give S C (T, O). 

The geodesies of S°(T, O) are described in the oriented graph (T, O) as fol- 
lows; let C and 1Z denote the matrices 

£=l I M R fl 



o l y V 1 1 

A closed path V of length A; on the graph may be described by starting at the 
midpoint of an edge, and then giving a sequence (iui, . . . , Wk), where each Wi 
is either I or r, signifying a left or right turn at the upcoming vertex. We then 
consider the matrix 

M v {k)=W 1 ...W k , (3) 

where Wj = C if w 3 ■, = I and Wj — TZ if Wj = r. 

The closed path V on T is then homotopic to a closed geodesic 7("P) on 
^(r, O) whose length length(7(P)) is given by 

2cosh( length ( 7(7?)) ) = tr(M P ). (4) 

We have to check what happen to the geodesies on S°(T, O) as we close the 
cusps and get 5 , °(r, O). Given a simple geodesic 7 on S°(T, O) let 7 be the 
imageof7inS c (r,£>). 

First 7 might be homotopicaly trivial. This can happen in two ways. One is 
that the orientation on the path is uniform (all "Left" or all "Right") and in this 
case the path is circling a cups and the length formula © gives since there is 
no geodesic in its homotopy class in S , °(r, O). The second possibility is that 7 
is nontrivial on S°(T, O) but bounds a disk in S' c (r, O), in this case the corre- 
sponding cycle on the graph disconnects the graph. 

Second, the image of two non equivalent geodesies on S C (T, O) might be- 
come homotopicaly equivalent on S C (T, O). The treatment of this case is similar 
to the second case above since, in this case the images of the two geodesies bound 
a cylinder in S' c (r, O) and therefore the cycles on the graph disconnect the graph. 

The probability that one or two cycles will disconnect the graph tends to by 
the following 
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Lemma 3.1. Let V be a cycle or union of two cycles on T n with length(V) < 
C log n then 

Prob„[V disconnect T n ] — > 0. 

The proof of this lemma is a straight forward application of the following 
results 



Theorem 3.1 (0). 



2 

Prob n [h{T) > —) —* 1 as n —* 00 • 



Therefore if "P disconnects T it induces a subgraph H C T with i>(i?) < 
Const * v(V) but it is known that 

Theorem 3.2 (|22l|6|). Let H be a graph such that v(H) < e(H). Then the 
expected number of copies ofH in a random cubic graph with n vertices is 0(n~ l ) 

Given a geodesic 7 on S , °(r, O). Let 7 be the image of 7 in S C (T, O) then if 
7 is not homotopicly trivial after closing all cusps, then there is a unique geodesic 
7 homotypical to 7 in S C (T, O). Under the large cusps condition we can compare 
the length of 7 and 7 



Theorem 3.3 (0). 

length{^) 
1 + e 



< length{^) < length^) 



Before computing the expected length of the geodesic corresponding to a cy- 
cle of length I on the graph we must deal with the question of how far the Bollobas 
estimate of Poisson distributions holds. It is clear that for very long cycles the dis- 
tribution is different lfl2l . The distribution of Hamiltonian cycles is known lfl2l . 
At the other extreme, for very short cycles the Poisson distribution is a very good 
estimate. Recently McKay, Wormald, Wysocka examine this problem [18| and 
gave a very detailed estimate for a more general problem. We will modify part 
of their results to suit our needs. In particular, their result covers the distribution 
of non-intersecting cycles. We will allow a very limited intersection, namely we 
allow for two cycles to touch once. In this case we can easily "pull" the corre- 
sponding path apart on the surface and thus the resulting geodesic on S C (T, O) 
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Figure 3: We cannot allow two intersections. 



will be non-itersecting. Note that allowing more then one component in the inter- 
section of two cycles will result in counting some geodesies more then once, as 
seen in|3]. 

We use the following to determine how long cycles may be so that the prob- 
ability that the intersection of any two cycles will have two or more components 
goes to as n — > oo. 

Theorem 3.4. Let K be a cubic graph, and C\> be the collection of all cycles of 
length r in K, with 3 < r < a. Let C C C\> x Cj, s < r, be all pairs of cycles 
(Ci, C2) with C\ f] C*2 7^ and C\ ^ G% In addition, suppose that the number of 
components p in G\ fl C 2 is at least 2. Given a cubic graph T with n vertices, 

E p(c\\jc 2 cv)<o m y, ^"(f (5) 

(Ci,c 2 )ec i>i,P>2 u >' v J 

where j is the number of edges in G\ fl Ci- 

Proof. Follows directly from formula 2.7 in ( |[T8lO . □ 

Summing for all pairs r,s < a, and assuming that a is chosen so that a 3 = 
o(n), we find that the probability that any cycle of length less than a has intersec- 
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tion with more than one component with a different cycle also of length less than 
a is 

/ 2 2«-lX 

O —H (6) 



n 2 



If we choose a so that 

2 2 "- 1 = o{n 2 ) 
then this probability will go to zero as n — > oo. Thus, 

2 20 - 1 = o(n 2 ) (7) 

2 2a-l = ^2-e (g) 

2a -1 = (2-e)log 2 n (9) 
a = (l-e')log 2 n (10) 

4 Lengths of Geodesies 

We consider cycles of length N in the graph, and the geodesies on the surface as- 
sociated with each cycle. The length of the geodesic associated to any given cycle 
can be computed using ©. Clearly these lengths will depend upon the orientation 
of the vertices, and for cycles of length N there are 2 N possible orientations. We 
need to understand the distribution of the lengths of the associated geodesies. 

We begin by considering the expected value and the distribution for the traces 
of the matrices M v from ©, where the paths V are some random set of paths 
of length N. We will do this by considering the individual matrix entries. Let 
M v {%) = WiW 2 ...Wi with i < N, and M v (0) = Id, the 2 x 2 identity matrix. If 



M v (i) 



a b 
c d 



then Mp(z + 1) is either 



or 



a + b b 
c + d c 

a a + b 
c c + d 

Consider the top row of M v {i). At step i = it is (10), and thus 1 and are the 
only values possible for either entry. At step i = 1, the top row is either (10) or 
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(11), and as i continues to increase, it is clear that the possible values are values 
from the previous step, or sums of neighboring values from the previous step. 
Thus, we can build a table of possible values for matrix elements by starting with 
1 and 0, putting the sum 1 in between, then putting the sums 2 and 1 in the new 
spaces, . . . 

1...2---1---1---0 (11) 

After i steps we have a sequence with 2 % + 1 entries. This sequence is one 
example of a Stern sequence ( ll20lfT3lfT4lfT5ll \ Stern sequences have many nice 
properties, and the moments at any step i can be computed directly. For example, 
each term, except for the initial terms a and b, shows up as part of two new terms 
in the next sequence, so if S(i) is the sum of terms at step i, then S(i + 1) = 
3S(i) — (a + b). Stern sequences, in the guise of the Stern-Brocot tree also show 
up in the work of Viswanath ([21 1) on random Fibonacci sequences, where he 
considered random products of the matrices 

, /0 1\ D / 1 
A= \l 1 J B= \l -1 

We need only consider values of diagonal elements of M v {i). At any given 
step, one diagonal element will change and the other will not, so it is clear that 
we need to understand both the distribution of the diagonal elements and the de- 
pendence between them to determine the distribution of values of the trace. The 
rows of the matrix associated with a path V will be generated by Stern sequences 
starting with (1 0) and (01), which we denote S u and Si. . These are simply reflec- 
tions of each other, and we omit the leading (or trailing) zero from the sequence. 
Counting from i = 0, we see that the sum of elements is 

i-l 

S(<) = J3 3< + 2 . 

5=1 

while the number of elements N(i) = 2\ hence the expected value at step i is 

We can compute the variance of by determining the sum of the squares of the 
elements in the sequence, S 2 (i), as follows (" |[T9lO . Consider three neighboring 
terms in the sequence 
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Setting 



and 



A(i) = X>i 

B(i) = J2xjX j+1 

one gets the recurrence relations 

A(i + 1) = 3A(i) + 2B{i) - 2 (13) 
B(i + 1) = 2A(i) + 2B(t)-2 (14) 
= 5A(i - 1) - 2A(i - 2) - 1 (15) 

A messy but straightforward computation gives 

2- I - 1 (l72 I (-5+v^) + (5+v^) I (-34+6v^) + (5-v / 17) I (-51+ll\/l7)) 



A(i) 



17(-5+vT7) 



using A[l] = 2, A[2] = 7. Similar calculations give the sums of nth powers of 
terms of the sequence. Thus we can compute any central moments of our Stern 
sequence, and the distribution of diagonal matrix entries. 

We must convert our knowledge of the moments of the diagonal entries into 
information about their sum. Since means are additive, we can compute the mean 
value of the trace using (fT2l . 

E _ 2 / E;=i3' + 2 \ _ 3 ' + l 

E Tr ,i-2\ - TI )-~^~ (16) 

Clearly a dependence exists between the value of the upper diagonal and lower 
diagonal elements x u and x\. Hence, determining the variance of the traces will 
require a computation of the covariance of the Stern Sequences for the upper and 
lower diagonal elements. The covariance can be computed using 

Cov(S h S u ) =< s t s u > - < Si>< S u > (17) 

where <> indicates the expected value. Equation (fT2l gives us the second term, 
so we need to compute the mean of products of associated elements in the two 
sequences, in other words the product of the diagonal entries for each matrix Mp 
corresponding to paths of length i. 
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We compute the product of diagonal matrix elements using the diagram in Fig- 
ureHJ which follows directly from an investigation of the matrices. After an initial 
choice of either C or 1Z, the diagonal product is 1 . Each number in a given row 
will generate two children, corresponding to a choice of C or 1Z. We determine 
the value of the children using two simple rules. If the path to the child is a contin- 
uation of a previous path from higher in the diagram, we add the same value that 
was previously added, otherwise we add the parents value to its' siblings, subtract 
one, and that becomes the addend for the new direction. The outer edges of the 
tree have addends of zero, which correspond to paths containing only C and only 

n. 



1 

+0/ +<i+i)-i=fi 

1 2 2 

+tV \ +(1+2)-1=+2/ \ +1 +1/ \ +{1+2)-1=+2 y v +0 

7\ / \ 7\ /\ 



Figure 4: Computing the products of corresponding elements of the Stern se- 
quence. 

As in the computation for the sums of squares, we get a recurrence relation 

C{i + l) = 5C(i) + 2C{t)-2 i - 1 (18) 

where 

= ^ Xu,i x l,i 
i=l 

andC(l) = 2, C(2) = 6. Thus, 

C{x) = ^2--- 2 (l72 2 -+ 1 -(-17+^)(5+v / ^) a: + (5-v / T7) a; (l7+v / T7)) (19) 

Substituting equations (fT2l and (fT9l into (fTTt gives 

Cova^ance = ^±g!!!!zg(g!h^ (20 ) 
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Looking at a plot of (127)1) in Figure |5J we see that while for short cycles there is a 
negative correlation between the diagonal elements, this quickly changes over to 
a positive one. In fact, we can compute the correlation between the two diagonal 
matrix entries, and see that 

lim Cor = *-C=- ~ .61 

n^oo 34 _ 6^17 




Figure 5: The covariance for the two diagonal elements of a matrix Mp(n). 



Knowledge of the covariance for the upper and lower diagonal elements of the 
matrices allows us to compute the variance for the traces, which gives us a basic 
understanding of the distribution of traces for cycles of length N. 

Variance = 4~ x (-1+2* -23* +a x -9* +{b-Vyf) x +(5+v/i7) x ) (21) 

Unfortunately, the lengths of geodesies on the surface corresponding to the given 
cycles is determined using ©, 

Tr f Tr -I- \/Tr 2 — 1 \ 
Length^) = 2 cosfT 1 (^-) = 2 In I —— — j . (22) 

When one has two independent random variables related by some function Y = 
f(X), the standard technique used to deduce information about the distribution 
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of the Y's given information about the distribution of the X's is to use a Taylor 
approximation, and evaluate at the expected value of X, E(X). 



Y 



f(c) + f(c)(X-c) + 



f{c){X-cf 



— - Letc = E(X) 
f'{E{X)){X-E{X)f 



(23) 



2 



Y 



f(E(X)) + f'(c)(X-E(X)) + 



2 



(24) 



One can now approximate E(Y) using (|24|> . The second term disappears, and 
since in the case we are considering /(X) rs ln(x), the third term gives roughly 



which approaches — oo exponentially. Thus, we get no lower bound even for the 
mean of the lengths as n increases. 

The following Lemma, which follows directly from the arithmetic geometric 
mean inequality, shows that the means of the lengths cannot be growing any faster 
than linearly in n. 

Lemma 4.1. For a positive random variable X 



We need to give a lower bound for the growth of the lengths. The following 
proposition shows that the lower bound is also linear. 

Theorem 4.1. The proportion of the standard Stern sequence which is growing 
exponentially in N as N — > oo approaches 1. 

Proof. First, recall that any pair of neighbors in the Nth Stern sequence represents 
a row in Aip(N). Next, note that in JAp(l) one row must be (1, 1), without loss 
of generality we will assume that this is the bottom row. Thus, at the expense of 
one step, we can consider the growth of the standard Stern sequence beginning 
with 1 1, which contains 2 k pairs after k steps. 

We may now consider a Bernoulli process defined as follows: for any pair 
in the kth Stern sequence, consider the four pairs which are its' children in the 
k + 2nd Stern sequence. These come from the four choices of turns possible, 
1Z1Z, TZC, CTZ, and CC Since we are considering the trace of Aip, we need 
only consider the right element in each pair. If our initial pair is a b, we must 



Var(X) 

E(xy 



E{\og(X)) < hg(EX) 
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a b 



:-; a+b d 



a 2ano a+b a 1 2b b 

Figure 6: Three out of four entries at least double every two steps 

consider two possibilities, a > b or a < b. The possibility that a = b will only 
occur in the last pair, which has vanishing probability. 

In Figure ©, we see that the four possible values for the diagonal element are 
2a + b, a + b, a + 2b, and b. If a < b, then one of the four is double the original 
diagonal element b, while if a > b, three of the four are double the original. 
As N — > oo, the probability that a < b and a > b approach 50%. Thus, defining 
success for our Bernoulli process as a doubling, the probability of success is 1/2 — 
e(N), with limjv^oc e = 0. 

Considering paths of length iV = 2k + 1, the Bernoulli process will have k 
steps. As N — > oo, the mean will approach N/2 and the probability that there will 
be at least k successes will approach one. Therefore the probability that any entry 
in the nth Stern sequence is at least 2 k approaches 1 as N — > oo. □ 

Keeping better track of the minimal growth using a table like Figure fallows 
us to give a numerical estimate for the lower bound on the growth rate. For exam- 
ple, if b < a, then of four choices for the second term in a pair of Stern sequence 
entries, two have at least tripled (2a + b and b + 2a). One then has a multinomial 
process, and if the number of steps corresponding to a single trial is allowed to 
increase, the lower bound also increases. 

Allowing five steps, one calculates that a lower bound for the growth factor of 
almost every element is given by 

237/32031/16^3/80^/40^1/80^1/160 ^ j 35592 

Such calculations quickly get too long to do by hand, and the return on investment 
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becomes minimal since the mean is growing like 1.5 . A simple Mathematica 
routine can push the calculations significantly farther, however; for example if the 
basic step size is 15, then the growth factor is approximately 1.43925. 
Combining this with d22l and Lemma |4~T1 we get 

Theorem 4.2. As N — > oo, the expected value of the length of the simple closed 
geodesies associated to cycles of length N, E N {1), is bounded by 

(logl.43925)iV < E N (l) < (log 1.5) JV. 

5 The Length of the Shortest Geodesic 

In this section we will use the results from previous section to estimate the length 
of the shortest closed geodesic on S C (T, O). The length of the shortest closed 
geodesic is an important geometric invariant since it is twice the injectivity radius. 
Let syst(S ,<7 ) (T, O) be length of the shortest closed geodesic on S C (T, O), we will 
show that : 

Theorem 5.1. 

2.809 < syst(S c )(r,0) < 3.085 

It is interesting to compare this result in with the the behavior of the largest 
embedded ball on(S c ) (T, O) which gives a linear growth to the largest embedded 
ball. 

We will start with the upper bound. As we have seen on a random graph the 
distribution of short cycles is independent on the size of the graph and it is Poisson 
distribution with mean A = £r where k is the length of the cycle. Therefore the 

probability of not having a k cycle is e - ^ and the probability of having at least 

one A;-cycle on the graph is (1 — ). 

As we transition from the graph T to the compact surface S C (T, O), by Lemma 
13. II there are only 2 orientations that will make the curve that corresponds to a cy- 
cle null homotopic. Hence the probability that there is a A;-cycle that gives rise to 

non trivial geodesic on S (T, O) isp(fc) = 2 2 k-2 (1 — e ~ ™ ) , and the probability 
that there is no A;-cycle that produces a geodesic on S C (T, O) will be (1 — p(k)). 
The expected value for the trace of a k cycle is tr{k)^^- and the expected value 
for the length of the geodesic is bounded from above by by 2 cosh _1 (^^). 

E(syst(S c )(r,0))<^(^)n( 1 -^)))2cosh- 1 (^ r ) 

k=3 j=2 
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= t - «-*) n (> - - «-*))><»*-0 

It is easy to see that this series converges rapidly and we can get a numerical 
estimate that the value is ~ 3.085 

To get a lower bound we can replace log E(tr(M)) by E(2 cosrT^^)) and 
use the rapid decay for the probability that a graph has large girth. We calculate 

k=2 ^ j=2 ~ ' 

and get the estimate of 2.809. 
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